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NOMENCLATURE 


a  Crack  depth 

A  Autofrettage 

c  Surface  crack  length 

C  Coefficient  in  Paris'  crack  growth  law 

E  Modulus  of  elasticity 

G  E/2(l+v) 

He  Correction  factor  defined  in  equation  (11) 

K  Stress  intensity  factor 

Kp  Opening  node  stress  intensity  factor 

Kpi  Sliding  node  stress  intensity  factor 

K  Non-dimens ionaiizing  stress  intensity  factor 

Kq  Non-dimensionalizing  stress  intensity  factor 

AK  Stress  intensity  factor  range 
Ji-n  Natural  logarithm 

m  Exponent  in  Paris'  crack  growth  law 
max  Maximum 

min  Minimum 

N  Number  of  loading  cycles 

p  Pressure 

Pil  Pure  bending 

PT  Pure  tension 

r  Radius 

R^  Tube  inner  radius 


iv 


R.2  Tube  outer  radius 

Ra  Tube  autofrettage  radius 

Rp  Tube  reversed  yielding  radius 

t  Theoretical 

Y  Yield  stress 

z  Complex  variable,  x  4-  iy 

a  Proportion  of  tension 

3  Proportion  of  bending 

Y  Tube  opening  angle  (radians) 

t  Parameter  plane 

0  Angular  coordinate 

k  3-4v  (plane  strain),  (3-v)/(l+v)  (plane  stress) 

X  Yield  stress  in  tension/yield  stress  in  compression 

v  Poisson's  ratio 

5  Parameter  plane 

o  Direct  stress 

t  Shear  stress 

<(>  Complex  stress  function 

t|>  Complex  stress  function 
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INTRODUCTION 


Fatigue  crack  growth  arising  from  the  cyclic  pressurization  of  thick-wall 
cylinders  tends  to  produce  radial  fatigue  cracks  emanating  from  the  bore.  A 
knowledge  of  the  crack  tip  stress  intensity  factor,  K,  is  necessary  in  order 
to  predict  the  fatigue  crack  growth  rate,  critical  length  and  lifetime  of  such 
cracks.  It  is  common  practice  to  produce  an  advantageous  stress  distribution 
by  autofrettage  (overstrain)  of  the  cylinder  in  order  to  alow  or  prevent  crack 
growth.  This  autofrettage  process  may  Involve  plastic  strain  throughout  the 
wall  thickness,  or  any  lesser  proportion  of  the  wall  thickness,  depending  upon 
the  degree  of  overstrain  applied  to  the  cylinder  by  over  pressure  or  by  an 
oversized  mandrel-swage  process. 

An  accurate  stress  intensity  solution  for  pressurized,  autof rettaged 
thick  cylinders  is  a  fundamental  requirement  for  crack  grcwth  rate  and  life 
prediction.  Some  of  the  work  relating  to  two  and  three  dimensional  K  solu¬ 
tions  is  reviewed  by  Tan  and  Fenner, ^  Of  the  various  types  of  solutions,  the 
errors  associated  with  collocation  and  integral  equation  solutions  are  of  the 

order  of  1  percent,  while  4  percent  would  be  more  typical  of  finite  element 

* 

and  boundary  element  methods.  However,  there  may  be  more  significant  uncer¬ 
tainties  in  crack  growth  and  life  prediction.  Some  key  factors,  and  the  fac¬ 
tors  considered  here,  are  the  shape  of  the  crack  during  the  fatigue  lifetime 
of  the  component,  uncertainty  over  the  exact  proportion  of  overstrain  in  the 
tube,  and  a  residual  stress  distribution  which  generally  does  not  conform  with 

*Tan,  C.  L.  and  Fenner,  R.  T. ,  "Stress  Intensity  Factors  for  Semi-Elliptical 
Surface  Cracks  in  the  Pressurized  Cylinders  Using  the  Boundary  Integral 
Equation  Method,"  Int.  J.  Fracture,  Vol.  16,  No.  3,  1980,  pp.  233-245. 


the  predictions  of  an  idealized  elastic-plastic  analysis,  particularly  the 
assumption  of  the  same  magnitude  of  yield  strength  in  tension  and  compression. 

In  this  report  it  is  proposed  that  accurate  two  dimensional  K  solutions 
may  be  modified  in  order  to  predict  stress  intensity  factors,  and  hence  crack, 
growth  rates,  for  semi-elliptical  cracks  in  pressurized  thick-wall  cylinders 
with  residual  stress  distributions.  Such  predictions  may  be  compared  with 
experimental  crack  growth  data,  and  some  measure  of  the  extent  of  overstrain 
may  be  obtained  by  cutting  autof rettaged  tubes  radially,  and  measuring  the 
angle  of  opening. 

CALCULATION  OF  STRESS  INTENSITY  FACTORS 

Stress  intensity  factors  for  the  plane  (two  dimensional)  geometry 
illustrated  in  Figure  1(a)  were  obtained  by  the  modified  mapping  collocation 
(MMC)  method.  This  method  is  described  in  detail  in  Reference  2.  Briefly, 
complex  variable  methods,  due  to  Muskhelishvili3  are  utilized.  Stresses  and 
displacements  within  a  body  are  given  in  terms  of  the  complex  stress  functions 
<t>(z)  and  i|j(z)  by: 

ax  +  Oy  »  4  Re  {4>’(z)}  s  (1) 

oy  -  ax  +  2i  Txy  -  2[zf’(z)  +  ♦•(*)]  (2) 


^Andrasic.,  C.  P.  and  Parker,  A.  P.  ,  "Weight  Functions  for  Cracked  Curved 
Beams,"  Numerical  Methods  in  Fracture  Mechanics,  D.  R.  J.  Owen  and  A.  R. 
Luxmoore  (Eds . ) ,  Proceedings  Second  International  Conference,  Swansea,  1980, 
pp.  67-82. 

^Muskhelishvili ,  N.  I.,  Some  Basic  Problems  of  the  Mathematical  Theory  of 
Elasticity,  Noordhoff,  1973. 


2G(u+iv)  -  k  <p(z)  -  zip  *  (z )  -  *(z)  (3) 

where  the  complex  variable  z  ■  x  +  iy,  x  and  y  are  the  physical  coordinates, 
primes  denote  differentiation,  and  bars  represent  the  complex  conjugate. 

Also: 

E 

G  - - 

2(  1+v) 

3-v 

k  -  3  -  4v  (plane  strain),  <  ■  -  (plane  stress) 

1+v 

where  G  is  the  shear  modulus,  E  the  elastic  modulus,  and  v  is  Poisson’s  ratio, 
while  the  resultant  force  over  an  arc  s  is: 

fl  +  if2  -  i  /  (Xn+iYn)ds  -  ij>( 2)  +  z  <j>'(z)  +  <Kz)  (4) 

s 

where  Xn  ds  and  Yn  ds  are  the  horizontal  and  vertical  components  of  force 
acting  on  ds. 

The  solution  of  the  cracked,  autof rettaged  cylinder  was  carried  out  on 
similar  lines  to  that  of  Tracy. 4  A  complex  mapping  transforms  straight  lines 
parallel  to  the  real  axis  in  the  5-plane  to  curved  lines  in  the  physical  (z) 
plane,  in  particular  the  real  axis  in  the  5  plane  is  mapped  to  an  arc  of 
radius  ,  centered  at  the  origin  in  the  z-plane,  Figure  1.  A  further  mapping 
is  introduced  which  maps  the  unit  semi-circle  plus  its  exterior  in  the  c  plane 
to  the  crack  plus  its  exterior  in  the  5-plane,  see  Figure  l(b,c).  The 
analytic  continuation  arguments  of  Muskhelishvili  are  used  to  ensure  traction- 


^Tracy,  P.  G. ,  "Elastic  Analysis  of  Radial  "racks  Emanating  From  the  Outer  and 
Inner  Surfaces  of  a  Circular  Ring,"  Engng.  Frac.  Mech. ,  Vol.  11,  1979,  pp. 
291-300, 


free  conditions  along  F'L'  and  J'G’  and  hence  FL  and  JG  in  the  physical  plane, 
Figure  1. 

For  certain  geometries  it  is  necessary  to  use  partitioning^  to  obtain  the 
desired  accuracy.  The  partitioning  of  the  cylinder  is  shown  in  Figure  l.  in 
general  each  region  has  its  own  complex  stress  and  mapping  function.  Since 
there  is  symmetry  about  the  imaginary  axis,  only  region  I  and  that  pare  of 
region  II  to  the  right  of  the  imaginary  axis  need  be  considered.  When 
partitioning  is  used  it  is  necessary  to  "stitch"  along  common  boundaries,  by 
imposing  equilibrium  and  compatibility  of  displacements. 

In  the  MMC  method  the  infinite  series  representations  of  the  stress 
functions  are  truncated  to  a  finite  number  of  terms.  Force  conditions  are 
imposed  at  selected  boundary  points  which  give  conditions  on  the  unknown 
coefficients  in  the  stress  functions.  Thus  ecch  boundary  point  produces  two 
rows  in  the  main  matrix  A,  and  two  corresponding  elements  In  the  boundary 
conditions  vector  >,  where: 

A  x  =  b 

and  x  is  the  vector  of  unknown  coefficients.  The  common  ("stitched")  boundary 
points  are  used  to  obtain  conditions  relating  the  unknown  coefficients.  Each 
common  boundary  point  gives  four  rows  in  A  and  four  corresponding  zeros  in  h. 
In  general,  A  is  a  matrix  of  i  rows  and  m  columns,  where  l  and  m  depend  upon 
the  number  of  boundary  points  and  unknown  coefficients  respectively.  It.  was 
found  that  convergence  is  generally  better  when  2m  <  t  <  2.6  m,  and  this 

^  And  rasic,  C.  P.  and  Parker,  A.  P. ,  "Weight  Functions  For  Cracked  Curved 
Beams,"  Numerical  Methods  in  Fracture  Mechanics,  D.  R.  J.  Owen  and  A.  R. 
Luxmoore  (Eds . )  ,  Proceedings  Second  International  Conference,  Swansea,  1980, 
pp.  67-82. 


compares  well  with  other  workers.’’  A  least-square  error  minimization 
procedure  was  used  to  solve  the  overdetermined  set  of  linear  equations. 

Knowing  the  coefficients  for  the  stress  function  In  the  cracked  region,  the 
crack  tip  stress  Intensity  factor,  K,  may  be  determined  from:0 

K  -  Kx  -  IKli  -  2(2m)1/2  Limit  (z-Zc)1/2  •  <t>'(z)  (5) 

z  *  zc 

where  Kx  and  Kn  are  opening  and  sliding  mode  stress  Intensity  factors 
respectively,  and  zc  is  the  location  of  the  crack  tip. 

We  now  consider  the  various  loadings,  and  combinations  of  loadings  which 
will  be  necessary  for  the  solution  of  the  pressurized,  autof rettaged  tube. 

LoadingA:  Internal  pressure  acting  in  bore  and  cracks. 

An  accurate  MMC  opening  mode  stress  Intensity  solution,  Kx ,  is  available 
for  this  configuration2  and  is  shown  as  the  upper  curve  in  Figure  2.  This  MMC 
solution  wa3  obtained  by  superposition  of  an  all-round  tension  field  on  the 
outer  boundary.  The  results  are  presented  in  dimensionless  form  as  Kx/Kp, 


where : 


2R22 

Kp  -  - n - 7-  p(Tta)1/2 

T  22—r1  2 ) 


and  R2  the  outer  radius,  Rq  the  inner  radius,  a  Is  the  crack  depth  and  p  is 
the  internal  pressure  acting  in  both  b>re  and  crack.  Also  shown,  for  the 


t’Fas<.n,  K.  1),,  "A  Review  of  Least  Squares  Methods  for  Solving  Partial 
differential  Equations,"  Int.  J.  Num.  Meth.  in  Engr.g.,  Vol.  10,  197b,  pp. 
1021-1046. 

^Sih,  G.  C. ,  Paris,  P.  C.  ,  and  Erdc.gan,  F.  ,  "Crack  Tip  Stress  Intensity 
Factors  for  Plane  Extension  and  Plate  Bending  Problems,"  J.  Appl .  Mech.,  Vol. 
29,  1962,  pp.  306-312. 

2  Bowie,  0.  L.  and  Freese,  C.  E. ,  "Elastic  Analysis  For  a  Radial  Crack  in  a 
ClrcuLar  Ring,"  Engineering  Fracture  Mechanics ,  Vol.  4,  No.  2,  1972,  pp. 
315-321. 
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purposes  of  comparison,  are  points  from  a  solution  due  to  Grandt^  based  on  an 
approximate  weight  function  method  which  was  also  employed  in  the  derivation 
of  additional  solutions  referred  to  in  the  next  section. 

Loading  B:  Ideal  autofrettage  residual  stresses  in  uncracked  tube. 

An  MMC  program,  based  on  the  formulation  outlined  in  this  section,  was 
used  to  calculate  the  stress  intensity  results  for  full  autofrettage  (100 
percent  overstrain),  shown  as  the  lower  curve  In  Figure  2.  This  was  based  on 
the  ideal,  elastic-plastic  solution,'^  where  the  distribution  of  hoop  stress 
(ay)  in  the  uncracked  tube  is  given  by; 

•  R^  R22 

o0  =  -Y  Jin  (R2/Ri)  [1  +  "7 - 7  (1  +  -r-)l  +  Y[  1  +  in  (r/Ri )  |  (7) 

R  2  2  ”  R  i 2  r2 

and  Y  Is  the  uniaxial  yield  strength  of  the  material  (Tresca's  criterion)  or 
1.15  x  yield  stength  (vor  'ses'  criterion).  Points  obtained  from  a  soLution 
by  Grandt  are  shown  for  comparison  purposes.  The  approximate  results  of 
Grandt"’  for  a  cracked  tube  of  the  same  dimensions  subjected  to  staady-state 
thermal  loading  were  modified  in  accordance  with  Reference  11  to  make  possible 
a  comparison  with  the  calculated  autofrettage  results.  The  modification  is 
based  on  the  fact  that  the  residual  stress  distribution  f of  100  percent  over- 


"lirandt  ,  A.  K.,  "Stress  Lntensity  Factors  For  Cracked  Holes  and  Rings  Loaded 
With  Polynomial  Crack  Face  Pressure  OLstr  ibutions  ,'*  Lnt.  .J.  Fracture,  Vo  1  . 
14,  1078,  pp.  R221-R229.  " 

^llill,  R.  ,  Hie  Mathemat  1  cal  Theory  of  Plasticity,  Clarendon  Press,  Oxford, 

19  50.  .  ““ 

'()Grandt,  A.  F.,  "Two  Dimensional  Stress  Intensity  Factor  Solutions  For 
Radially  Cracked  Rings, "  AFML-TR-7 5-1 2 1 ,  Air  Force  Materials  Lab,  Wright 
Patterson  AFB,  Ohio,  1975. 

‘'Parker,  A.  P.  and  Farrow,  .J.  R.  ,  "Stress  Intensity  Factors  for  Multiple 
Radial  Cracks  Emanating  Form  the  Bore  of  an  Autof rettaged  or  Thermally 
Stressed  Thick  Cylinder,"  Er.gng,  Frac.  Mech.  ,  Vol.  14,  1981,  pp.  237-241. 
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strain  Is  Identical  to  that  for  steady  state  thermal  loading,  apart  from  a 
simple  multiplying  constant.  Agreement  is  generally  within  5  percent. 

A  particularly  important  feature  of  the  MMC  formulation  outlined  at  the 
beginning  of  this  section  is  the  accuracy  of  the  results  at  shallow  crack 
depths.  Many  available  results  are  not  quoted  for  a/CKo-^!)  <  0.1,  or  even 
0.2,  and  are  therefore  of  limited  use  for  life  prediction  purposes.  The  Lnset 
in  Figure  2  shows  the  calculated  results  for  100  percent  overstrain  In  the 
range  0  <  a/(R2~k[)  0.1,  and  indicates  good  convergence  to  the  limiting 

value  of  1.12.  Since  as  much  as  80  percent  of  gun  tube  lifetime  may  be 
expended  in  this  range  it  is  clearly  very  important  to  seek  accurate  results 
at  shallow  crack,  depths.  Ft  is  apparent  that  this  solution  provides  the 
required  accuracy.  The  results  for  full  autofrettage  are  presented  in 
dimensionless  form  as  Kj/ KA  where: 

2Rz2 

KA  »  [1  -  2n  f  K 2 / K i )  ( — -) |  Y  ( n a ) 1 / 2  (8) 

R?.  Ri 

Loading  C:  Internal  pressure  and  ideal  autofrettage. 

The  total  stress  intensity  Ln  a  pressurized,  fully  autof rettaged  (100 
percent  overstrain)  tube  is  given  h  the  superposition  of  the  results  In 
Figure  2,  thus: 

Kful  1  uuiolrott  f  pressure  =  Rp  *V\  1 > 

where  Kp  is  the  stress  intensity  with  pressure  in  bore  and  crack,  and  KA  is 
t. Vie  stress  Intensity  contribution  due  to  the  100  percent  overstrain  residual 
stress  field  acting  alone. 
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In  the  event  that  the  tube  has  been  subjected  to  less  than  100  percent 
overstrain,  the  plastic  flow  during  the  autofrettage  process  will  extend  to  a 
radius  Ra,  and  the  stress  intensity  factor  in  this  case  is  given  by: 

Y  Y 

^partial  autofrettage  ”  [  1  +  -  £h(R2/Ra)  -  -  ( R22~Ra2  ) / 2R22  I  +  > 

1  p  P 

+  pressure 

a  <  Ra  -  Ri  ( 10) 

where  Ra  can  be  obtained  from  the  expression:^ 

Y  9  9 

p  -  YHoCRa/R!)  +  - -  (R22~Ra2) 

2R  2 

Results  for  0,  30,  60,  and  100  percent  overstrain,  with  Y/p  -  3.53,  R2 / H 1 
=  2.0,  based  on  the  superposition  of  results  given  in  Figure  2  are  shown  in 
Figure  3.  These  curves  indicate  the  very  significant  reduction  in  stress 
intensity  as  a  result  of  the  ideal  autofrettage  process.  Indeed,  for  the 
particular  value  Y/p  -  3.55,  100  percent  overstrain  causes  a  negative  total 
stress  intensity  (i.e.,  crack  closure)  for  crack  depths  up  to  a/(R2-Ri)  * 
0.08. 

CORRECTION  FOR  THE  CRACK  SHAPE  EFFECTS 

Thus  far  we  have  ignored  the  effect  of  crack  shape  on  stress  intensity, 
assuming  a  through  crack.  We  now  consider  the  crack  to  be  semi-elliptical, 
semi-major  axis  c,  depth  a,  Figure  4.  In  order  to  modify  the  two-dimensional 
results  to  account  for  crack  shape  we  employ  two  sets  of  work,  namely  the 
extensive  results  for  semi-elliptical,  cracks  in  a  flat  plate  under  tension  or 

^Parker,  A.  P,  and  Farrow,  J.  R. ,  "Stress  Intensity  Factors  for  Multiple 
Radial  Cracks  Emanating  From  the  Bore  of  an  Autof rettaged  or  Thermally 
Stressed  Thick  Cylinder,"  Engng.  Frac.  Mec.h,  ,  Vol.  14,  1981,  pp.  237-241  . 
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beading,  due  to  Newman  and  Raju,^  and  the  limited  results  for  a  semi- 
elliptical  crack  in  a  pressurized  thick  cylinder,  due  to  Tan  and  Fenner^  and 
Atluri  and  Kathlresan.^ 

First,  consider  the  flat  plate  containing  a  semi-elliptical  crack.  From 
the  results  given  in  Reference  13,  we  may  define  a  correction  factor  for  the 
plate  in  pure  tension  Hp^  given  by: 

Vr 

•Vr  -  : — 

KpT 

where  Kp-p  is  the  stress  intensity  factor  for  semi-elliptical  crack  in  a  plate 
in  tension  and  Kpx  is  the  stress  intensity  factor  solution  for  a  straight- 
fronted  through  crack  given  in  Reference  15.  Curves  of  the  correction  factor 
Hp'f  are  shown  in  Figure  5(a). 

Similar  correction  factors  for  the  flat  plate  under  pure  bending  Hpg  are 
defined  by: 

*pb 

!1pb  *  r— 

where  Kpg  is  the  stress  Intensity  factor  for  a  semi-elliptical  crack  In  a 

""  \ 

plate  in  pure  bending  and  Kpg  is  the  solution  for  a  straight-fronted  through 
crack.  ^5  Curves  of  the  correction  factor  Hpjj  are  shown  in  Figure  3(b). 


^Tan,  C.  L.  and  Fenner,  R.  T.  ,  "Stress  intensity  Factors  for  Semi-Elliptical. 
Surface  Cracks  in  Pressurized  Cylinders  Using  the  Boundary  integral  Equation 
Method,"  Int.  J.  Fracture,  Vol.  16,  No.  3,  1980,  pp.  233-245. 

^Newman,  J.  C.  and  Raju,  t.  S.,  "Analyses  of  Surface  Cracks  in  Finite  Plates 
Under  Tension  or  bending  Loads,"  NASA  TP  1578,  1979. 
l^Atulrl,  N.  and  Kathi  resan,  K. ,  "3D  Analyses  of  Surface  Flaws  in  Thick- 
Walled  Reactor  Pressure-Vessels  Using  Displacement-Hybrid  Finite  Element 
Method,"  Nuclear  Engng .  and  Design,  Vol.  51,  1979,  pp.  163-176. 
l^Rooke,  D.  P.  and  Cartwright,  D.  J. ,  Compedium  of  Stress  Intensity  Factors, 
Her  Majesty's  Stationery  Office,  London,  1976. 
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It  is  proposed  that,  at  shalLow  crack  depths,  the  correction  factors 
applicable  to  the  thick  cylinder  may  be  obtained  by  appropriate  superpositions 


of  Hpx  and  Hpu,  given  by. 


H(j  -  aHp-p  +  (SRpu  (11) 

where  the  multiplying  factors  a  and  8  are  obtained  by  calculating  the 
proportions  of  tension  and  bending  in.  the  uncracked  tube  which  act  over  the 
prospective  crack  line.  In  order  to  test  this  hypothesis,  consider  the  case 
of  a  tube  subjected  to  internal  pressure.  In  this  case  the  crack-line  loading 
comprises  hoop  stresses  due  to  the  internal  pressure^2  plus  a  contribution 
from  the  pressure,  p,  which  infiltrates  the  crack,  thus  the  total  crack  line 
loading  is: 


00 


P(  1  + 


nr 


R22-Rj2 


*22 

( 1  +  -r-) ) 


as  shown  in  Figure  6(a).  Using  a  straight  line  approximation,  for  the  range 
0  <  a/(R2-Ri)  <  0.2,  the  proportion  of  tension  loading  is  given  by  a  ■ 
1.72/2.66  while  the  proportion  of  pure  bending,  0  -  0.94/2.66,  ’hence  the 
correction  factor  1^  is  determined  as 

He  *  .645  Hpx  +  .355  Hpjj 

Correction  factors  determined  on  this  basis  are  plotted  in  Figure  6(b) 
for  the  case  a/c  ■  0.8.  Also  shown  are  the  equivalent  correction  factors 
obtained  for  the  same  configuration  by  Tan  and  Fenner, ^  using  boundary 


Tan,  C.  L,  and  Fenner,  R,  T. ,  "Stress  Intensity  Factors  for  Semi-Elliptical 
Surface  Cracks  in  Pressurised  Cylinders  Using  the  Boundary  Integral  Equation 
Method,"  Int.  .J.  Fracture,  Vol.  16,  No.  3,  1980,  pp.  233-245. 

12Parker,  A.  P. ,  "Stress  Intensity  and  Fatigue  Crack  Growth  in  Multiply- 
Cracked,  Pressurized,  Partially  Autof rettaged  Thick  Cylinders,"  Fatigue 
of  Engng .  Materials  and  Structures  (in  press).  ~ 
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element  inethods.  Agreement  Is  within  2  percent  at  a/(R2~Ri)  ■  0.2. 
deeper  cracks  the  difference  between  plate  and  cylinder  results  becomes 
appreciable  because  of  the  differences  in  constraint,  discussed  in  the 
upcoming  paragraph.  It  is  important  to  emphasize  that  the  proportion  of  the 
tube  in  which  the  vast  majority  of  component  lifetime  is  expended  is  welL 
represented  by  the  model,  i.e.,  for  a/CRy-Ri)  0.3.  The  equivalent 
correction  factors  for  K  at  the  point  4>  -  0  (the  free  surface  of  the  plate) 
are  shown  in  Figure  6(c)  and  exhibit  similar  characteristics  to  those  of 
Figure  6(b). 

The  significant  difference  between  the  cylinder  and  plate  shape  factors 
in  Figure  6(b)  is  attributed  to  the  significantly  greater  constraint  of  a 
cylinder  compared  with  a  plate,  particularly  for  straight-fronted,  deep 
cracks.  This  difference  can  be  demonstrated  by  comparing  the  K  for  a 
straight-fronted,  a/(R2-Rl)  ■  0.8  crack  in  a  pressurized  cylinder  with  R2/R1 
2.0  with  the  K  for  a  straight-fronted,  a/(R2“Ri)  “  0.8  crack  in  a  plate  with 
approximately  the  same  combination  of  tension  and  bending  loading  as  that  of 
the  pressurized  cylinder.  The  K  for  the  cylinder  can  be  obtained  from 
Reference  7  or  from  Figure  2  and  is  shown  in  Table  I.  The  K  for  the  plate 
wLth  cylinder  loading  Is  obtained  using  the  expressions  for  pure  tension  and 
pure  bending  as  follows: 

Opiate  “  kRT  +  K-PB 

^  Bowie,  0.  L.  and  Freese,  C.  E.,  "Elastic  Analysis  For  a  Radial  Crack  In  a 
Circular  Ring,"  Engineering  Fracture  Mechanics,  Vol.  4,  No.  2,  1972,  pp. 
315-321.  —  — 

^Rooke,  D.  P.  and  Cartwright,  D.  J.,  Compedlum  of  Stress  Intensity  Factors, 
Her  Majesty's  Stationery  Office,  London,  1976. 
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The  tension  and  bending  stresses  in  Kpx  and  Kpg  are  determined  as  opT  =  2.0  p 
aiid  opg  -  0.5  p,  by  using  a  linear  approximation  of  the  entire  cylinder 
loading  plot  shown  in  Figure  6(a).  The  K  for  the  plate  with  a  straight  crack 
is  several  times  that  of  the  cylinder,  which  indicates  that  the  plate  is  much 
less  constrained.  Therefore,  the  increase  in  constraint  corresponding  to  a 
semi-elliptical  rather  than  a  straight  crack  will  be  much  larger  for  a  plate 
than  for  a  cylinder.  This  results  in  the  much  smaller  for  a  plate  than  for 
a  cylinder,  as  observed  in  Figure  6(b).  When  the  final  comparison  is  made  in 
Table  l,  K  for  deep  semi-elliptical  cracks  in  a  cylinder  and  a  similarly 
loaded  plate  are  about  the  same,  as  might  be  expected. 

TABLE  I.  COMPARISON  OF  K  VALUES  FOR  DEEPLY  CRACKED  CYLINDER  AND'  PLATE 


K  for  a/(R2-Ri)  - 

r—  -------------  ■  -  -  -  ■  ■ 

^  from  Fig.  6(b) 

K  for  a/(R2~Rl)  ■ 

0.8 

for  a/(R2~Ri)  - 

0.8,  a/c  -  0.8 

straight  crack 

0.8,  a/c  -  0.8 

semi-elliptical 

crack 

pressurized 

cylinder 

6.99  p(a)!/2 

0.42 

2.9  p(a)l/2 

plate  with 
approx, 
cylinder 
loading 

46. 7  p(a)l/2 

0.07 

3.3  p(a)l/2 

In  order  to  obtain  a  range  of  correction  factors  applicable  to  thick 
cylinders,  the  procedure  outlined  by  Eq.  (11),  Figure  6,  anu  related 


discussion  was  used  to  obtain  correction  factors  from  Newman's  work  for  the 
case  of  internal  pressure,  and  full  autofrettage  (100  percent  overstrain). 
These  results  for  relatively  short  cracks  were  combined,  using  engineering 
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judgment,  with  those  of  Tan  and  Fenner  for  longer  cracks  to  obtain  correction 
factors  for  use  over  a  wide  range  of  crack  length.  The  results  for  pressure, 
Lp,  are  shown  in  Figure  7(a)  and  7(b),  for  <|>  »  n/2  and  0.0  respectively,  and 
for  autofrettage,  La,  are  shown  in  Figure  7(c)  and  7(d)  for  <j>  *  it/ 2  and  0.0 
respectively.  In  the  case  of  pressure,  Figure  7(a)  and  7(b),  we  may  compare 
in  each  figure  with  the  solutions  due  to  Atluri  and  Kathiresan^  at  two 
points,  namely  a/c  ■  1.0,  a/t  »  0.5  and  0.8.  For  $  “  tt / 2  agreement  is  within 
8  percent,  however,  for  <j>  »  0  the  correction  factors,  differ  by  25  percent, 
reflecting  wide  disagreement  between  workers  on  K  solutions  at  the  free- 
surface. 13,14 

It  is  thus  possible  to  calculate  the  stress  intensity  factor  for  any 
combination  of  pressure,  ideal  partial  autofrettage  and  crack  shape  by 
correcting  the  stress  intensity  factors  presented  in  Figure  2  in  accordance 
with  Eq,  (10)  and  correction  factors  Lp  and  La  to  obtain: 

Y  Y 

^partial  autofrettage  “  ^p  t  ^  -  &h(R2/Ra)  “  “  ( R2?‘-Ra2 )/ 2R22  1  Kp 

+  pressure  (3D)  p  p 


+  lA  Ka 


(12) 


^Newman,  J.  C.  and  Raju,  I.  S.,  "Analyses  of  Surface  Cracks  in  Finite  Plates 
Under  Tension  or  Bending  Loads,"  NASA  TP  1578,  1979. 
l^Atulri,  N.  and  Kathiresan,  K. ,  "3D  Analyses  of  Surface  Flaws  in  Thick-Walled 
Reactor  Pressure-Vessels  Using  Displacement-Hybrid  Finite  Element  Method," 
Nuclear  Engng.  and  Design,  Vol.  51,  1979,  pp.  163-176. 
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EFFECTS  OF  NONIDEAL  RESIDUAL  STRESS  DISTRIBUTION 


The  residual  stress  distribution  predicted  by  Eq.  (7)  may  not  occur  in 
practice. 16,17  Non-ideal  Bauschinger  effects  and  uncertainty  over  the  exact 
amount  of  autofrettage  on  unloading  may  produce  a  different  stress  distribu¬ 
tion.  In  earlier  v*>rkl®  a  constant  reduction  factor  of  0.7  was  applied  to  the 
autofrettage  contribution  to  stress  intensity  factor  to  account  for  primarily 
the  Bauschinger  effect. 

In  the  present  vrerk  it  is  proposed  that  these  effects  may  be  considered 
in  a  different  manner.  If  an  unflawed  autof rettaged  tube  is  cut  along  a 
radius,  it  will  spring  apart,  the  theoretical  angle  of  opening  yt  being  given 
for  100  percent  overstrain  and  using  the  von  Mises'  yield  criterion  as:l® 


•100% 

Yt 


8irY 

/lE 


(13) 


In  the  case  of  partial  overstrain,  the  total  moment  acting  over  the  cut 
ends  is  reduced  by  a  factor  F.  Details  of  the  calculation  of  F  are  contained 
in  an  Appendix.  The  theoretical  angle  of  opening  for  partial  overstrain 
is  given  by: 

-  100% 

Yt*  "  pYt  04) 


l^Underwood,  J.  H.  and  Throop,  J.  F.,  "Residual  Stress  Effects  on  Fatigue 
Cracking  of  Pressurized  Cylinders  and  Notched  Bending  Specimens."  presented 
at  Fourth  SESA  International  Congress,  Boston,  MA,  May  1980. 
l^Milligan,  R.  V.,  too,  W.  H. ,  Davidson,  T.  E. ,  "The  Bauschinger  Effect  in  a 
High-Strength  Steel,"  Trans.  ASME,  J.  Bas.  Engineering,  Vol.  88,  1966,  pp. 
480-488. 

IBparker,  A.  P.  and  Farrow,  J.  R. ,  "On  the  Equivalence  of  Axisymmetric 
Bending,  Thermal,  and  Autofrettage  Residual  Stress  Fields,"  J.  Strain 
Analysis,  Vol.  15,  No.  1,  1980,  pp.  51-52. 


A  graphical  representation  of  F  for  a  tube  having  R2/R1  "  2.0  is  shown  in 
Figure  8.  The  angle  (and  moment)  ratio  F  between  partial  and  100  percent 
overstrain  does  not  vary  much  with  R2/R1.  In  fact,  for  tubes  in  the  range  1.8 
<  R2/R1  <  2.2  the  deviation  from  the  curve  in  Figure  8  is  only  1  percent,  and 
of  course  it  still  goes  asymptotically  to  the  limit  of  1.0  at  100  percent 
overstrain,  and  to  zero  at  0  percent  overstrain. 

We  propose  that  a  comparison  of  the  ratio  of  measured  opening  angle  to 

100% 

ft  with  the  ratio  F  provides  an  indication  of  the  non-ideality  of  the 
residual  stress  distribution  in  actual  autofrettage  cylinders.  The  data 
points  shown  in  Figure  8  make  this  comparison.  Each  point  represents  an  autn- 
frettaged  steel  cylinder  of  the  type  described  in  Reference  16.  A  section  wan 
removed  from  each  cylinder,  and  the  opening  angle  was  measured  using  an 
optical  comparator.  The  ratio  of  measured  angle  to  100  percent  overstrain 
theoretical  angle  was  calculated  using  the  measured  value  of  yield  strength,  Y 
for  each  tube. 

The  Important  features  of  the  comparison  between  experiment  and  theory  in 

Figure  8  are  the  following,  (a)  The  experimental  results  are  generally  near 

% 

or  above  the  theoretical  curve  for  relatively  low  overstrain  and  generally 
below  the  curve  for  high  overstrain'.  This  can  be  explained  by  the  Bauschinger 
effect  since,  for  high  overstrain  and  the  associated  large  amount  of  tensile 
yielding,  the  Bauschinger  effect  would  result  in  significant  reverse  yielding 
and  less  than  expected  residual  stress  and  opening  angle,  (b)  The  two  sets  of 

^Underwood,  J.  H.  and  Throop,  J.  F. ,  "Residual  Stress  Effects  on  Fatigue 
Cracking  of  Pressurized  Cylinders  and  Notched  Bending  Specimens,"  presented 
at  Fourth  SESA  International  Congress,  Boston,  MA,  May  1980. 
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experimental  results  designated  by  dashed  lines  indicate  less  than  expected 
opening  angle  and  residual  stress  with  increasing  R2/R1.  This  also  is 
consistent  with  the  Bauschinger  effect,  using  a  similar  rationale  as  with  (a), 
that  is,  larger  R2/R1  result  in  more  tensile  yielding,  more  reverse  yielding, 
and  less  than  expected  residual  stress. 

Also  shown  in  Figure  8  are  the  measured-to-theoretical  angle  ratios  from 
the  30  percent  and  60  percent  overstrained  tubes  of  the  experimental  work 
described  here.  Note  that  the  average  of  the  two  angle  measurements  for  60 
percent  overstrain  is  about  0.6  of  the  theoretical  value.  This  would  suggest 
that,  at  least  as  a  first  approximation,  the  contribution  of  overstrain  to  the 

total  K  could  be  reduced  by  this  0.6  factor,  and  a  shorter  than  expected 

\ 

fatigue  life  would  result.  This  is  discussed  further  in  the  upcoming 

comparison  of  experimental  and  theoretical  results. 

No  attempt  is  made  to  relate  the  two  sets  of  experimental  results  in 

Figure  8,  that  is  the  results  designated  by  a  dashed  line  for  50  percent  to 

100  percent  overstrained  tubes  and  the  results  from  the  30  percent  and  60 

percent  overstrained  tubes.  The  reason  for  this  is  that  the  methods  used  for 

% 

overstraining  were  quite  different  for  the  two  sets  of  results.  A  mandrel 
swaging  process  was  used  in  the  former  case  and  a  hydraulic  overpressure 
process  was  used  in  the  latter  case.  The  different  overstraining  methods 
could  result  in  different  amounts  of  Bauschinger  effect  in  the  two  sets  of 
experimental  results. 


LIFE  CALCULATIONS 


The  fatigue  growth  rate  of  cracks  subjected  to  cyclic  loading  may  be 
expressed  in  terms  of  Paris'  law: ^ 


—  -  C(  AlOm  (15) 

dN 


where  da/dN  is  the  fatigue  crack  growth  per  loading  cycle,  C  and  m  re 
empirical  constants  and  AK  is  the  ringe  of  stress  intensity  defined  by: 

*  l^aax  ”  ^min  (K®in  *  0) 

AK  *  K^ax  (^min  ^  0) 

and  Kqjqx  and  K^n  are  respectively  the  maximum  and  minimum  values  of  stress 
intensity  during  the  loading  cycle.  Note  that  the  possibility  of  "over¬ 
lapping"  or  touching  of  the  crack  surfaces  at  some  point  on  the  crack  line 
remote  from  the  crack  tip^O  is  not  considered  In  this  report.  During  the 
lifetime  of  a  particular  cracked  cylinder  the  crack  will  propagate  from  some 
initial  depth  ai  to  some  final  depth  af,  where  af  is  generally  the  total  wall 
thickness  of  the  cylinder,  R2-R1.  In  order  to  predict  the  fatigue  life, 
equation  (15)  is  rearranged  to  give: 


,af  da 

I  - -  Nf  -  Ni 

at  C( AK)n 


(16) 


^Paris,  p.  c.  and  Erdogan,  F.,  "A  Critical  Analysis  of  Crack  Propagation 
Laws,"  Trans  ASME,  J.  Bas.  Engng.,  Vol.  85,  1963,  pp.  528-534. 

20parker,  A.  P. ,  "Stress  Intensity  Factors,  Crac..  Profiles,  and  Fatigue  Crack 
Growth  Rates  in  Residual  Stress  Fields,"  presented  at  ASTM  Symposium  on 
Residual  Stress  Effects  in  Fatigue,  Phoenix,  AZ,  May  1981. 


Crack  growth  (a  versus  N)  predictions  are  made  for  two  examples  which  are 
near  the  extremes  of  crack  geometry  encountered  in  thick  cylinders,  a  single 
nearly  straight-fronted  crack  and  a  single  semi-circular  crack.  The 
predictions  are  compared  with  ultrasonic  crack  growth  measurements  from 
cylinders  in  which  internal  radius  is  90  mm,  R2/R1  is  2.0,  the  cyclic 
pressurization  is  zero  to  331  MPa.  The  cylinder  material  is  ASTM  A723  forged 
steel,  with  yield  strength  of  1175  MPA,  -40°C  Charpy  impact  energy  of  34  J, 
reduction  in  area  of  50  percent. 21  These  properties  can  be  slightly  different 
after  overstrain,  due  to  the  plastic  strain,  which  can  be  up  to  1  percent  at 
the  inner  radius  of  a  100  percent  overstrained  cylinder.  Considering  this 
small  amount  of  plastic  strain  relative  to  reduction  in  area,  no  significant 
effect  on  fatigue  life  is  expected  as  the  result  of  the  material  property 
changes  due  to  the  overstrain  process. 

In  general,  the  integral  of  Eq.  (16)  was  evaluated  numerically  using 
Simpson's  rule.  For  he  particular  steel  employed  in  the  experimental  crack 
growth  rate  wrk,  the  measured  constants  are: 

C  ■  6.52  x  10”^2  ,  m  ■  3.0 

for  crack  growth  in  meters  per  cycle  and  AK  in  MPa 
(a)  Single,  nearly  straight-fronted  crack 

Figure  9  shows  a  versus  N  predictions  based  on  the  K.  results  presented 
in  Figure  3  for  a  single,  straight-fronted  radial  crack,  initial  depth  6.4  ram. 
The  solid  lines  in  Figure  9  are  predictions  for  zero  and  30  percent  over- 


2^Throop,  J,  F,,  "Fatigue  Crack  Growth  in  Thick-Walled  Cylinders,"  Proceedings 
of  the  National  Conference  on  Fluid  Power,  Fluid  Power  Society,  Chicago, 

19T2 ;  pp.  115-T3T. 


18 


strain.  The  dashed  Line  shows  experimental  results,  originally  reported  in 
Reference  21,  for  a  tube  with  zero  nominal  overstrain  in  which  a  single  notch 
was  cut  using  electro-discharge  machining  to  a  depth  of  6.4  ram  and  a  half 
surface  length,  c,  of  254  mm. 

For  this  example  it  is  possible^  to  integrate  Eq .  (16)  directly  if  a 
simple,  shallow  crack  K  expression  is  used,  that  is,  one  with  Kp/ Kp  at  a 
constant  value.  From  Figure  2,  Kp/Kp  ■  1.05  ac  a/R2-Ri  ■  0.2  is  a  reasonable 
choice,  considering  that  most  of  the  cylinder  life  is  expended  at  relatively 
low  values  of  a.  Doing  so  and  combining  with  Eqs.  (6)  and  (12)  gives  the 
following  expression  for  fatigue  life  of  a  tube  with  no  overstrain: 


Nf  -  Ni 


1 

/af” 


Kp  2R22 

c(r-  (*)1/2  (-— o — 2)Lp  p]3 

K>,  R22-Ri2 


(17) 


For  values  of  a^,  C,  R2 ,  Rj  ,  and  p  in  this  example  and  setting  Lp  ■  1.0,  the 
result  is  the  dotted  line  in  Figure  9,  very  close  to  the  more  general 
analysis.  For  the  conditions  of  this  comparison,  the  simpler  analysis  of  Kq, 
(17),  although  less  rigorous  and  general,  is  adequate  and  easier  to  use. 

The  lack  of  agreement  between  the  two  analytical  predictions  and  the 
experiment  would  be  improved  if  the  shape  factor,  Lp  wa.3  significantly  less 


^Throop,  J.  F.,  "Fatigue  Crack  Growth  in  Thick-Walled  Cylinders,"  Proceedings 
of  the  National  Conference  cn  Fluid  Power,  Fluid  Power  Society,  Chicago, 
1972,  pp.  115-131. 

22underwood,  J.  H.  and  Throop,  J.  F.,  "Surface  Crack  K-Estimates  and  Fatigue 
Life  Calculations  in  Cannon  Tubes,”  Part-Through  Crack  Fatigue  Life 
Prediction,  ASTM  STP  687,  J.  B.  Chang,  Ed.,  American  Society  for  Testing  and 
Materials,  1479,  pp.  195-210. 
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than  unity.  The  Initial  shape  is  described  by  (a/c)i  m  0.025  and  the  final, 
shape  is  (a/c)f  ■  0.18.  Referring  to  Figure  7(a)  and  considering  that  most 
of  the  cylinder  life  is  expended  with  low  values  of  both  a/c  and  a/R2~Ri,  Lp 
is  estimated  to  be  between  0.95  and  1.0.  So  even  with  this  factor  to  the 
third  power  in  Eq.  (17),  it  cannot  account  for  the  differences  in  Figure  9. 

(h)  Single,  semi-circular  crack 

Figure  10  shows  a  versus  N  predictions,  based  on  the  K  values  presented 
in  Figure  2,  as  modified  in  accordance  with  Eq.  (12),  utilizing  correction 
factors  Lp  and  La  from  Figures  7(a)  and  7(c)  respectively,  for  a/c  »  1.0.  Tin* 
K  values  are  representative  of  the  stress  intensity  at  the  deepest  |x>lnt  of  a 
semi-circular  crack  of  initial  depth  6.4  mm  in  a  thick  cylinder,  and  therefore 
implicitly  assume  that  the  crack  retains  its  semi-circular  shape  during  the 
life  of  the  tube.  This  assumption  appears  to  be  justified  on  the  basis  of 
experimental  observations.^  The  predictions  are  made  for  zero,  30  percent, 

40  percent,  45  percent,  and  50  percent  overstrain,  while  the  dashed  lines 
indicate  experimental  results^  for  the  same  initial  crack  depth  for  actual 
tubes  with  zero,  30  percent,  and  60  percent  nominal  overstrain. 

For  this  example  it  is  also  possible  to  obtain  a  closed  fora  expression 
for  fatigue  life,  analogous  to  Eq ,  (17).  By  combining  Eqs.  (6),  (8),  and  (In; 
it  can  be  shown  that  the  ratio  of  the  K  expression  for  an  autofruttagi 1  and 
pressurized  cylinder  to  that  of  a  pressurized  cylinder  with  no  autofret  tag**  I  ; 
the  following: 

^Underwood,  J.  H.  and  Throop,  J.  F.,  "Residual  Stress  Effects  on  Fatigue 
Cracking  of  Pressurized  Cylinders  and  Notched  Bending  Specimens,"  presented 
at  Fourth  SESA  International  Congress,  Boston,  MA,  May  1980 
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(18) 


Kp+A  Y  r2  r1  rA2  -  <SRi 2  -  (1-6)R22 

Kp  1  +  p  ^  +  R2  +  2R22  1 


Ka/Ka 

whare  6  Is  defined  as  - s-  ,  the  ratio  of  the  K  solutions  in  Figure  2. 

Kp/Kp 


Using  Eq.  (18)  in  a  life  expression  in  the  form  of  Eq.  (17)  gives  Nf  -  Nj  * 


2R22  Y  r2  r1  rA2 

(fl)1/2(~:r - -)(1  +  -  [£n( — )  +  5*n  —  + - 

R22“R12  P  rA  r2 


( 1-6)R22  - 
- ])Lp]3 

(19) 


which  should  give  a  good  estimate  of  fatigue  life  for  an  autof rettaged , 
pressurized  cylinder  in  which  a  shallow  semi-elliptical  crack  dominates  the 
life.  Results  from  this  shallow  crack  expression  are  shown  in  Figure  10  for 
0  percent,  30  percent,  and  50  percent  overstrain.  A  mean  L  ■  0.53  was  used, 
since  Lp  and  are  close  to  this  value  for  a/R2-Ri  *  0.2,  see  Figures  7(a) 
and  7(c).  The  predictions  of  the  shallow  crack  analysis,  Eq.  (19),  are  close 
to  those  of  the  more  general  analysis.  The  lack  of  agreement  between  the  two 
analytical  predictions  and  the  experiments  is  discussed  in  .the  upcoming 
section. 


DISCUSSION  AND  CONCLUSIONS 

Life  prediction  for  cracked  tubes  requires  an  accurate  knowledge  of 
stress  intensity  factor  at  short  crack  depths.  In  this  report  two-dimensional 
solutions  for  cracked  tubes  with  various  amounts  of  residual  stress  (over¬ 
strain)  were  obtained  by  use  of  the  modified  mapping  collocation  technique, 
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which  gives  good  conveno-m-  ■  ...•  .  i  . . •  1  analytic  solution  at  very  small 

crack  depths. 

The  solutions  were  extended  i.»  include  three-dimensional  semi-elliptical 
cracks  in  thick  cylinders  by  <  ■<  t  -nl  o  i  ng  the  proportion  of  tension  and  bending 
in  the  cylinder  and  apply  In-,  op.unje  '.lie  correction  factors  determined  for 
cracked  plates  in  tension  no.  1  M  e.-,  Comparison  with  Limited  available 
three  dimensional,  solution-.  '  >•  •hi-  k  vlinders  Indicates  good  agreement  at 
short  to  medium  crock  I  n,-,  :  i.  . ,  "  :  *•  n.d,  wherein  most  of  the  fatigue 

lifetime  is  expended.  The  mm--  "  h  i  tupted  is  however,  not  restricted  to 
thick  cylinders,  and  m  iv  lee.  <  ,  :■  c,ipl  Icntion  in  fracture  Mechanics 

Design.  Indeed,  the  thick  ■  I  in'.  .1-1 y  he  a  sort,  of  limiting  configuration  in 
which  maximum  errors  wool  1  hi  ..n,  >•  i  noted  because  of  the  Large  degree  of 
restraint  offered  by  the  doubly  connected  body. 

The  calculations  of  llicil'ic  '  ->i  mm  -autof  rettaged ,  pressurized  tubes 
with  nearly  straight-fronted  •  -nl  -circular  cracks  predict  lives  of  about 
one  quarter  of  those  deiei  ••  e  1  ,  .  i  men  tally.  Possible  explanations  for 

this  effect  arc: 

a.  The  possibility  fhol  !o!  I  iit-essure  doe:-:  not  infiltrate  the  crack. 

This  is  not  considered  lit.  j-  >  :>e  expeiiment.il  work  was  deliberately 

designed  to  avoid  l  hi  ••  e:  •  "mild  h-  tor  uvuo  noticeable  in  the 

autof  rettage  case  ii  ii  w  ■,  •  h.-c  tuse  of  the  compressive  stress 

closing  the  crack. 

b.  Residual,  st.rcr.ses  In  1  ee  •i.in--outo>f rettage  tube.  Again,  this  is  not 

considered  likely.  When  .  m-.intoi  i  ettaged  tubes  do  not  exhibit 

any  tendency  Co  spring  op-  t. ,  •  iegl.igi.hle  gross  residual  stresses. 

77 


c.  Multiple,  small,  semi-elllptlcal  cracks  along  the  notch  boundary. 

This  is  considered  a  probable  explanation.  Such  a  form  of  multiple  crack 
growth  would  result  in  slower  overall  crack  growth  than  that  of  a  single 
crack,  until  the  individual  semi-elliptical  cracks  linked  to  form  a  single 
continuous  crack  front. 

Lifetimes  calculated  on  the  basis  of  the  'ideal'  residual  stress  fleLd 
indicate  extreme  sensitivity  to  the  amount  of  overstrain.  Note  in  Figure  10 
that  an  increase  from  30  to  50  percent  overstrain  increases  the  predicted  Life 
by  a  factor  of  five,  while  an  increase  from  30  to  60  percent  overstrain 
increases  measured  life  by  a  factor  of  1.4.  So  there  is  clearly  much  less 
increase  in  life  due  to  increased  overstrain  than  would  be  expected  from  the 
calculations.  We  believe  this  effect  is  directly  related  to  the  less  than 
expected  opening  angle  measured  from  cylinders  with  60  percent  overstrain, 
discussed  in  relation  to  Figure  8.  The  preferred  explanation  for  the 
deviations  of  opening  angle  and  lifetime  from  calculated  values  is  the 
reduction  in  the  residual  stress  field  due  to  reversed  yielding  near  the  inner 
radius  caused  by  the  hauschinger  effect.  This  would  manifest  itself  as  a 
reduction  in  opening  angle,  particularly  at  large  overstrains,  and  a 
relatively  larger  reduction  in  lifetime,  since  most  of  the  lifetime  is 
expended  at  shallow  depths  at  which  the  hoop  stresses  are  reduced  most 
significantly. 
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Figure  3.  Stress  intensity  factors  for  a  single,  straight-fronted  crack  in 
pressurized  thick  cylinder  with  0  percent,  30  percent,  f>0  percen 
and  100  percent  overstrain. 
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Figure  7.  (a)  Crack  shape  factors,  Lp,  for  pressure  in  bore  and  cracks, 

4>  -  ir/2;  two  data  points  from  Reference  14,  a/c  «  1.0. 

(b)  Crack  shape  factors,  Lp,  for  pressure  in  bore  and  cracks, 
4>  -  0 ;  two  data  points  from  Reference  14,  a/c  =  1.0. 

(c)  Crack  shape  factors,  La ,  for  100  percent  overstrain,  ■}>  - 

(d)  Crack  shape  factors,  La,  for  100  percent  overstrain,  4>  = 


Theoretical  Results: 
Eq.  14 
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Ratio  of  opening  angle  for  partial  overstrain  to  the  theoretical 
value  for  100  percent  overstrain. 
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Figure  9.  Fatigue  crack  growth  (a  versus  N)  in  thick  cylinder  with  single, 
straight-fronted  crack;  experimental  results  and  predictions. 


Figure  10. 
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Fatigue  crack  growth  (a  versus  N) 
semi-circular  crack;  experimental 


in  thick  cylinder  with  single 
results  and  predictions. 


APPENDIX 


OPENING  OF  CUT  TUBES  WITH  PARTIAL  AUTOFRETTAGE  -  MO  REVERSED  YIELDING 

The  residual  >tresses  in  an  autof  rettaged  r.ube ,  internal  radius  R^  , 
external  radius  R2,  autof rettage  radius  R/  -jre:^ 

R2  r2Z 

001  -  -P  +  Y(i  +  4n(r/Ri)  -p(— 5- - r  [1  +  ---])  ,  Ri  <  r  <  RA  (Al) 

R2  “R1  r 


The  total  moment  acting  over  any  r  i 4  r».  1  cut  is  given  by: 

rR2 

M  *>  j  00  •  r  •  dr 
R1 


rA  r2 

M  -  J  O0i  r  •  dr  +  J  002  r  •  dr  (A4) 

R1  rA 


Since  the  opening  angle  of  fully  autof rettaged  tubes  of  any  radius  ratio, 
r2/r1>  provided  reversed  yielding  does  not  occur,  is  given  by  8ttY//3e,^2  and 
the  amount  or  opening  is  proportional  to  the  applied  bending  moment,  it  is 


A^Hill,  R. ,  The  Mathematical  Theory  of  Plasticity,  Clarendon  Press,  Oxford, 
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possible  to  produce  a  non-dimensionalized  plot  of  theoretical  angle  of  opening 
for  a  cut  tube  with  partial  autof rettage ,  for  any  radius  ratio.  The  curve  for 
R2/R1  *»  2.0  is  shown  in  Figure  8  of  the  main  report.  Deviations  from  this 
curve  for  1.8  <  R2/R1  <  2.2  are  less  than  1  percent  of  the  maximum  opening. 


OPENING  OF  CUT  TUBES  WITH  PARTIAL  AUTOFRETTAGE  -  WITH  REVERSED  YIELDING 

If  the  material  of  the  tube  has  a  reduced  yield  strength  in  compression 
o’  -XY,  it  may  undergo  reversed  yielding  out  to  a  radius  Rq  after  removal  of 
tn«s  autof iettage  pressure,  p.  In  this  case  the  residual  stresses  ara:^-^ 

003  r,  -AY  (1  +  4n(r/Ri))  »  Rl  <  r  <  RD  ( AS] 


rd 


004  “  "fp  -  (1+X)  YJln  (--)}  —  r - r  [1  +  -;- 

Rl  R22“Rd2  r2 


-p  +  Y  (1  +'An(r/Ri>)  ,  RD  <  r  <  RA  (A6) 

YRA2  RD  rD2  r22 

095  -  [ - r  -  (p  -  (1+X)  Yin  ( — ) }  — r - -]  [1  +  ---]  ,  RA  <  r  <  R2  (A7) 

2R22  Rl  r22-Rd2  r2 


where  Rp  is  calculated  in  an  iterative  fashion  from: 

D-2_R^2 


p  R2  -RD 

-  -  (1+X)  [ - - —  +  2n( Rp/ Ri ) 

Y  2R22 


Once  again,  the  total  moment  acting  over  any  radial  cut  is  given  by 
equation  (A4) ,  thus  when  reversed  yielding  occurs : 

Rd  ra  r2 

M  -  /  003  r  •  dr  +  /  094  r  •  dr  +  /  005  r  .  dr  (A9) 

Rl  Rd  Ra 
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Tubes,  Some  Numerical  Methods  and  Results,"  US  Army  Numerical  Analysis  and 
Computers  Conference,  Huntsville,  AL  (1981)  (in  press). 


The  reduction  in  opening  arising  from  a  Bauschinger  effect  equivalent  to 
X  ■  0.5  can  be  calculated,  leading  to  a  plot  similar  to  Figure  8  of  the  main 
report.  The  maximum  reduction  in  opening  is  approximately  8  percent  as  a 
result  of  this  magnitude  of  Bauschinger  effect.  Deviations  from  this  figure 
do  not  exceed  1  percent  for  1.8  <  R2/ R-l  <  2.2. 
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